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Abstract. Recent studies have indicated that the coarse grained dynamics of a large 
class of traffic models and driven-diffusive systems may be described by urn models. 
We consider a class of one-dimensional urn models whereby particles hop from an 
urn to its nearest neighbor by a rate which decays with the occupation number k of 
the departure site as (1 + b/k). In addition a diffusion process takes place, whereby 
all particles in an urn may hop to an adjacent one with some rate a. Condensation 
transition which may take place in this model is studied and the (b, a) phase diagram 
is calculated within the mean field approximation and by numerical simulations. A 
driven-diffusive model whose coarse grained dynamics corresponds to this urn model 
is considered. 



PACS numbers: 02.50.Ey, 64.75.+g, 89.40.+k 
1. Introduction 

Ordering and condensation transitions in one- dimensional systems far from thermal 
equilibrium have been studied extensively in recent years [H 121 El- It has been repeatedly 
demonstrated that unlike systems in thermal equilibrium, driven one-dimensional 
systems whose dynamics does not obey detailed-balance can be ordered even when 
the dynamics is local and noisy. 

By making a correspondence between phase separation in one-dimensional systems 
and condensation in urn models, two mechanisms which allow for phase separation in 
driven systems have been suggested. Urn models are simple lattice models defined on a 
ring geometry, where each site can either be vacant or occupied by one or more particles. 
The first mechanism is described in terms of the Zero Range Process (ZRP) In this 

model particles hop between nearest neighbor lattice sites with rates Uk which depend 
only on the number of particles k at the departure site. If the rates decay to zero in 
the large k limit, or if the rates decay slowly enough to a finite value, a condensation 
transition takes place, whereby a single lattice site becomes macroscopically occupied as 
the density is increased beyond a critical value. It has been suggested that the coarse- 
grained dynamics in a broad class of one-dimensional driven models can be described 
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by a ZRP with rates which at large k decay as = uo^ (1 + b/k) E]- I n this case 
phase separation takes place only if b > 2. 

A second mechanism is described in terms of the of the Chipping Model (CM) 
[3 |H]. The dynamics of this model involves two processes: chipping, where a single 
particle hops to a nearest neighbor site at a constant rate u ; and diffusion, where all 
particles in a site hop together to an adjacent site with rate a. This model can be 
viewed as a ZRP with constant hopping rate (b = 0), extended to include diffusion 
processes as well. Mean-field analysis of this model indicates that this model exhibits 
a condensation transition at a critical density. This result remains valid beyond the 
mean-field approximation as long as the chipping process is symmetric [Zj. It has also 
been shown that if the chipping process is biased either to the left or to the right, no 
condensation transition takes place 8\. 

Recently it has been suggested that the coarse-grained dynamics of certain cellular- 
automata traffic models can be modeled by the asymmetric CM j^j. Within this picture, 
the coarse-grained evolution of traffic models is described in terms of domain dynamics, 
which essentially involves asymmetric chipping and diffusion processes. It has thus been 
concluded that no phase separation transition should be expected in this class of models, 
and that jamming phenomena take place as a broad crossover process rather than via 
a sharp phase transition. On the other hand, there exist models which may be related 
to traffic, in which the rate of chipping a particle from a domain is not constant, but 
rather depends on the domain size. Such an example is the Bus Route model ^H] > where 
an approximate description assigns hopping rates to buses that decay as a function of 
the distance to the next bus ahead. It is thus of interest to combine the features of the 
chipping model and the ZRP processes and consider urn modes which exhibit chipping 
and diffusion processes with occupation dependent chipping rates. 

In this paper we consider in details a class of urn models which incorporate both 
diffusion and chipping processes with occupation (k) dependent chipping rates of the 
form = 1 + b/k. In Section 2 we introduce the generalized ZRP and analyze 
its behavior both within mean-field approximation and by numerical simulations. In 
Section 3 we consider a simple driven diffusive model, and show that its coarse grained 
dynamics is described by this generalized ZRP. Conclusions and summary are presented 
in Section 4. 

2. Zero-Range Processes with diffusion 

The generalized zero-range process is defined on a lattice of M sites, with periodic 
boundary conditions, occupied by N = <pM particles. The dynamics is defined through 
the rates by which two nearest neighbor sites containing k and m particles, respectively, 
exchange particles: 

diffusion: (k, m) + m, 0) 

chipping: (k, m) ^+(k + 1, m — 1) , (k, m) 1 — ^ h (k — 1, m + 1) . 
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where the rate Uk depends on the number of particles in the departure site, 

b 

u k = l + -. (1) 

Thus the model is specified by three dynamical parameters, a, b and q, and by the 
average occupancy <fi. The limit a = corresponds to a ZRP with no diffusion; the 
limit 6 = recovers the chipping model. The parameter q controls the spatial bias, with 
q = 1/2 being the symmetric case. 

2.1. Mean-Field analysis 

We first analyze the model within mean-field approximation, where correlations between 
sites are neglected. Let pk be the probability of a given site to be occupied by k particles. 
The mean-field evolution equations for pk are given by 

= a(l -p ) 2 + ^iPi(l -Po) ~ (A-Wipi)po ( 2 ) 
dp n ~ 1 

-^j- = -2ap n (l- p ) + a^2p k p n „ k + [A(p n _i - Pn) + U n +lPn+l ~ ^nPn] , 

fe=l 

with n > 1 and A = J2T=i u kPk- 

To solve these equations in the steady state, where all time derivatives vanish, let us 
define a generating function g(s) = Y^=iVnS n /n. Multiplying the nth equation in (J2J) 
by s n and summing over n > 1 one obtains 

as 2 g\sf - 2 as(l - Po )g'(s) + (1 — s)(l — A s)g'(s) (3) 

+ b^—^g{s) - Ap (l -s) + a{l - p ) 2 = . 
s 

It is assumed that at criticality p k has the asymptotic large k behavior pk ~ l/k T . 
One is interested in evaluating r as a function of a and b. For the model to exhibit 
a condensation transition one needs r > 2. Within the mean- field approximation the 
model exhibits condensation transition already at b = 0. We therefore expect such 
transition to take place for any b > (an assertion which has been checked numerically), 
and hence assume in the following that r > 2. This form of pk implies that g(s) is well- 
defined for \s\ < 1, and is singular at s — 1. Taking s = e~ e , with e > 0, we find that 
for non-integer r the singular part of g(s) has the following behavior for e < 1, 

/oo e -en 
^TT rfn ~ eT ' ( 4 ) 

where for the purpose of extracting the leading singularity we replace the sum by an 
integral. For integer r the singularity of g(s) is of the form e T loge. In the following 
analysis we assume for convenience that r is non-integer. The results derived in this 
analysis apply to the integer r case as well. We thus make the anzats 

g(s) = (oq + axe + a 2 e 2 + ...)+ e T (6 + he + b 2 e 2 + ...)+..., (5) 

where the terms in the first brackets represent the regular part of g(s), while the terms 
in the second brackets correspond to the singular part, with r being its leading power. 
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Figure 1. The exponent r as a function of b for a = 0.2. Solid line is given by the 
solution to the mean- field equations. Boxes are results of numerical simulations of the 
mean- field dynamics with M = 1000, = 3. An estimation for t is obtained by fitting 
Pk to k~ T over the range 10 < k < 100. An example for such a fit in the case b = 3 is 
given in the inset. 



Clearly non-leading singularities with powers r' > t can exist in this expansion as well. 
In the following we proceed by inserting this anzats into (jUJ), and systematically solving 
this equation order by order in e. From the equations for terms of order 0, 1,2 and r 
one readily derives expressions for both the critical occupancy <p c and po in terms of A : 

, A (1-A) 2 

Note that non-leading singular terms cannot contribute to terms of these orders, and 
they thus need not be considered. 

In order to determine A and r one needs to distinguish between two cases : 

(i ) r > 3- - One proceeds by considering the orders 3,4 and r + 1 terms in ©• 
Since r > 3, non-leading singularities do not contribute to these terms as well. A 
straightforward calculation yields 

A = 2a + t+1 _^ (?) 



r = 9b 



36- ^(47 + 362 + 3-66 -4(2a + 6+l) v ^) 1/2 *, (8) 



where 7 = 4a (a + b + 1) + {b — 2) 2 . Inserting (JJJ) into (0) one gets the critical 
occupancy <p c in terms of the model parameters. This solution holds for values of b 
larger than 6*, defined by r(b*) = 3. Note that in the limit b — > 00 one has r — > b . 
This is a result of the fact that at large values of b diffusion becomes rare and the 
ZRP result is approached. 

{ii) 2 < r < 3 — In this interval we find that for b < b* (for which r < 3) r becomes 
independent of b, with r = 5/2. To analyze this case we consider the 2r — 2 term 
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Figure 2. Mean- field critical occupancy 4> c as a function of 6 for several values of a. 
Solid lines are given by © and J7J for r > 3. Dots are obtained numerically from 
© for r < 3, as explained in the text. The a — line is given by the ZRP solution, 
(f> c = 1/(6 — 2). The <0> symbols indicate the value of 6* on each curve. At 6 = the 
critical density is given by the CM solution <fi c = •v/l + 1/a — 1. The value of 6* as a 
function of a, as obtained in mean-field, is given in the inset. 

in (JHJ). This term is given by aT 2 b\e 2T ~ 2 , which does not vanish since by definition 
b 7^ 0. In this r interval non-leading singular terms in g(s) cannot contribute to 
this term. To see this, let r' > r correspond to a such sub-leading term. Its leading 
contribution to (j3J) is of the order e r +1 , which is of higher order than e 2r ~ 2 . This 
analysis implies that 2r — 2 has to be an integer. Within this r interval, the only 
possibility is r = 5/2. 

The parameter A, and thus the critical occupancy, cannot be determined 
perturbatively in this r interval. It is determined by the two boundary conditions 
imposed on Q, namely g(0) = and g(l) = (A — 1 +po)/b, with p given by (Pi- 
Integrating 01 numerically, one identifies the value of A which satisfies the boundary 
conditions. 

In summary, we find that within the mean- field approximation r = 5/2 for 
b < b*(a), while r is a continuous function of b for b > b*(a). The exponent r 
exhibits a discontinuity at b* (a) (see figure The critical occupancy C is given as a 
function of b in figure EJ For b > b*(a) the curve is deduced from and while for 
b < b*(a) it is obtained numerically by finding the value of A for which the solution of (J3J) 
satisfies the appropriate boundary conditions. One observes that although the exponent 
t exhibits a discontinuity at b*, the critical occupancy is a continuous monotonically 
decreasing function of b. Note that increasing either a or b results in decreasing the 
critical occupancy. This is due to the fact that condensation is favored by both the 
diffusion process and by the size-dependency of the chipping. The value of b* as a 
function of a is depicted in the inset of this figure. 
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Figure 3. The exponent r as a function of b, obtained from numerical simulations. 
Solid lines are r — 2 corresponding to asymmetric CM (6 = 0), and r = b corresponding 
to ZRP (a — 0). Simulations were performed on systems of size M — 1000 and 
occupancy <fi — 3. 

The r{b) curve found above corresponds to the mean-field approximation in the 
thermodynamic limit. It is of interest to examine how this behavior is manifested in 
finite systems. This could give a useful insight for the model with nearest-neighbor 
exchange, where no analytic results are available, and where one has to resort only to 
numerical simulations of finite systems. To this end we carried out numerical simulations 
of the model, with a modified dynamics for which the mean-field equations yield the 
correct steady state. This modified dynamics is defined by the same rates as those of the 
original model, except that the two sites exchanging particles in each dynamical move 
are chosen at random, and are not necessarily nearest neighbors. The results obtained 
for t from a simulation of a system of size M = 1000 are presented in figure ^ It 
is observed that the qualitative behavior derived in the previous section is recovered. 
Namely, r assumes the value 5/2 for small b, while it continuously varies with b at 
large b. In the vicinity of the critical parameter b*, large crossover effects dominate 
the dynamics, and the numerical curve seems to deviate from the analytical one. This 
is a finite size effect, and one needs far larger systems in order to recover the true 
asymptotic exponent r in this region, as given by the exact solution found above. We 
carried out similar calculations for larger systems and found that indeed the trend is to 
move towards the analytical curve. 

2.2. Numerical simulations 

In the case of nearest-neighbor exchange no analytical solution for the steady-state 
distribution function is available, and one has to resort to numerical simulations. The 
exponent r(b) for the totally asymmetric case (q = 0), as extracted from studies of 
systems of size M = 1000, is given in figure H3 It is readily seen that the numerical 
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results are consistent with r = 2 at small b and with a continuously varying r at large 
b. This indicates that as in the CM (b = 0) the generalized model does not exhibit a 
condensation transition for small b > 0. Condensation is obtained only above a certain 
threshold of the parameter b. The precise nature of the curve in hard to determine 
numerically, due to the finite size effects which dominate the dynamics in the vicinity 
of b*, as for the mean- field dynamics. Thus, for example, accurate determination of b* 
is not possible, although it seems close, or equal, to 2. It is also not clear from the 
simulations whether or not r(b) is discontinues at b*. 

We also carried out numerical studies of the symmetric case, q = |. For b = it 
is known that r = 5/2, as in mean-field [Zj. We find that r remains locked at 5/2 for 
sufficiently small b, while r varies continuously above this value for large b. 



3. Corresponding one-dimensional driven diffusive model 

We now demonstrate that the generalized zero-range process considered in the previous 
section could be relevant to the coarse-grained dynamics of certain driven-diffusive 
models. To this end we introduce a simple driven-diffusive model, and analyze it using 
the generalized ZRP model discussed above. This model is related to a class of models 
whose coarse-grained dynamics has been described in terms of the ZRP jSJ|B]. In contrast 
to those models, the model presented here exhibits also domain diffusion. The model 
evolves in discrete time (namely through a parallel update scheme), and is therefore 
related to cellular automata models introduced to study traffic flow [TT| IT2"]. 

The model is defined on a lattice of L sites. Each site can either be vacant (0), or 
occupied by positive (+) or negative (— ) particle. The dynamical moves are carried out 
in parallel by two consecutive update steps. 

Step 1: 

Step 2: 

■0-0 
^000 -0 
^00 - 00 . 

Here x is a particle of either type and u + a < 1. The model allows for single and 
double site jumps, and it exhibits features which do not exist in models with only 
single-site jumps. Note that while all transition probabilities in step 2 are symmetric 
under particle-type exchange and left-right reflection, this symmetry is broken by the 
transitions of step 1. At high density the system has states in which all vacancies are 
isolated, and are thus deterministic. In this case ergodicity is broken. 

Consider now the coarse-grained dynamics of this model. Let us define a domain 
in this model as a sequence of particles, interrupted only by isolated vacancies. These 
vacancies move deterministically to the left. The evolution of a domain can be described 
in terms of two processes: (a) a chipping process, in which a positive (negative) particle 
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Figure 4. Space time configurations of the driven model starting from random initial 
configuration. A system of size L — 300 is occupied by 100 particles of each type, and 
the model parameters are a — 0.45, u = 0.5 and p = 0.1. Positive particles are colored 
black, negative particles are in gray, and vacancies are in white. 

leaves a domain of size k to the right (left) with rate w£ (w^)] and (b) a diffusion 
process, in which a vacancy penetrates the domain from its right (with probability a) 
and advances deterministically to its left, thus shifting its center of mass one unit to 
the right. In figure 0] a space-time configuration of the model is given. One readily 
observes the evolution of domains through chipping and diffusion, and the coalescence 
of domains upon contact. 

We now consider the current which flows through a domain. The size k of a domain 
is defined as the number of particles it contains, ignoring the vacancies in it. Since 
vacancies move deterministically through a domain, one may project the vacancies out of 
the domain dynamics, and consider only the dynamics of the charges. This dynamics is 
in fact described by the totally asymmetric exclusion process (TASEP) ^Hj- Typically, a 
domain is asymmetric, composed of an unequal + and — charge densities. In its maximal 
current state, obtained for large a and u, the current through such a domain takes the 
form j k ~ 1 + b/k with b = 3/2 [HI EH]. We thus conclude that w± = w± (1 + b/k), 
with u>+ 7^ w x some constants, and b = 3/2. 

The coarse-grained domain dynamics of this model is thus related to the generalized 
ZRP picture of Section |2l with q ^ | and b = 3/2. Since our numerical data (figure EJ) 
suggests that b* > 3/2, we expect the domain size distribution pk to take the form 
Pk ~ k~ T , with t = 2. Within the domain length accessible in numerical studies one 
could not get an accurate estimate for the exponent r. However, numerical simulations 
of the model indicate that r is larger than its zero- diffusion limit (3/2), and is close or 
equal to 2. Note that the correspondence between this model and the zero-range process 
can only be made in the region of the model parameters where the domains are in the 
maximal- current state, and where the absorbing states are not reached. 

Before concluding let us note that nearest-neighbor attractive interaction between 
particle of the same species could result in larger values of b, possibly reaching b > b* ~ 2. 
In this case the dynamics within a domain becomes identical to that of the KLS model 
[ini HI]. Here = (1 + b{rj)/k), where rj is the strength of the interaction For 
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large values of rj one expects the exponent r to increase beyond r = 2. This should be 
accompanied by a phase separation at high densities. 

4. Summary and Discussion 

In this work we generalize a class of zero-range processes with hopping rates of the 
form to (1 + b/k) to include diffusion. In particular the exponent r associated with the 
occupation distribution p k ~ k~ T , and its dependence on the parameter b, is studied. 
Both mean-field approximation and numerical analysis show that for sufficiently small b 
the exponent r is b independent, and assumes its b = value. In particular, in the mean- 
field approximation r = 5/2 while for asymmetric nearest-neighbor chipping numerical 
simulations indicate that r = 2. Above a critical value of b the exponent r becomes a 
continuously increasing function of b. 

This model could be relevant for describing the coarse-grained dynamics of certain 
driven diffusive models, and possibly traffic models. A particular example of a driven 
diffusive system for which this is the case is presented. 

The fact that r = 2 in the asymmetric chipping-model has been interpreted as an 
indication that jamming in traffic models do not take place via a sharp phase transition 

The analysis presented in this paper indicates that this result is rather robust, 
and could hold even for a more general class of traffic models, for which the chipping- 
like process exhibits a weak dependence on the length of the domain (b < b*). Such 
dependence may be of relevance for traffic models which allow overtaking. 
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